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ABSTRACT 


In  this  report  the  reviced  technique  for  designing  a  tactical  rocket 
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has  been  replaced  by  a  digital  flight  computer.  This  report  presents  a 
design  technique  for  designing  a  digital  attitude  control  system  for  the 
PERSHING  weapon  system.  The  design  goal  is  to  match  the  digital  system 
response  to  that  of  the  original  analogue  system. 
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1.  Introduction 


A  technique  under  development  by  Professor  Benjamin  C.  Kuo  of 
the  University  of  Illinois  was  applied  to  the  design  of  a  tactical  rocket 

attitude  control  system.  The  applicability  of  the  technique  has  been 
increased  and  the  analysis  used  ir«  Reference  1  has  been  extended 

accordingly 

In  this  report  the  revised  technique  nas  been  extended  and  applied 
to  the  redesign  of  a  continuous  control  system.  The  original  PERSHING 
\  weapon  system  attitude  control  system  was  implemented  by  an  analogue 

\  flight  computer.  This  computer  has  been  replaced  by  a  digital  flight 

{  computer.  This  report  presents  a  design  technique  for  designing  a  digi- 

5  tal  attitude  control  system  for  che  PERSHING  weapon  system.  The  design 

|  goal  is  to  match  the  digital  system  respon' p  .u  that  of  the  original 

I  analogue  system. 


2.  Design  Technique 

The  theoretical  basis  is  presented  in  References  1  and  2.  Only 
the  basic  design  technique  will  be  described  here. 

Assume  that  a  linear  continuous  system  (Figure  1)  is  described  by 
the  state  variable  equations, 

x(t)  =  A  x(t)  +  B  m(t)  (1) 

and 

m(t)  =  E(0)  r(t)  -  G(0)  x(t)  ,  (2) 

where  x(t)  is  a  vector  representing  the  state  of  the  system,  and  A,  B, 
G(0),  and  E(0)  represen  time -invariant  matrices.  The  design  goal  is  to 
replace  this  system  with  a  digital  version  (Figure  2)  whose  error  signal 


■'"Seltzer,  S,  M.,  Digital  Redesign  of  a  Rocket  Attitude  Control 
System,  U.  S.  Army  Missile  Command,  Redstone  Arsenal,  Alabama;  January 
1971,  Report  No.  RG-TN-71-1. 
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Seltzer,  S.  M. ,  Digital  Redesign  of  a  Rocket  Attitude  Control 
System.  Part  II,  U.S.  Army  Missile  Command,  Redstone  Arsenal,  Alabama, 
November  1971,  Report  No.  RG-TN-71-5. 
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m(t)  is  identical  to  the  continuous  system  error  m(t).  The  discrete  data 
system  im.de  1  that  is  used  to  replace  the  continuous  model  may  be  represented 
analytically  by 


x(t)  =  A  x(t)  +  B  m(t) 


and 


m(t)  =  E (T)  r(t)  -  G(T)  x(t) 


(3) 


(4) 


where  G(T)  and  E(T)  represent  time- invariant  matrices.  As  indicated  in 
References  1  and  2,  matrices  G(T)  and  E(T)  are  found  by  matching  corre¬ 
sponding  terms  of  truncated  Taylor's  series  expansions  of  x(t)  and 

x(t)  yielding: 


G(T)  «  G(0) 


+  T 


SG(T)] 

ST  1 

JT  =  0 


and 


where 


and 


E (T)  «  E (0) 


+  T 


SECT) 

ST 


iT  = 


> 


0 


Sg(t)~ 

ST 


T  =  0 


\  G(0)  [A  -  B  G (0) ] 


SE(T)~ 

ST 


T  =  0 


\  G(0)  BE (0) 


(5) 


(6) 


(7) 


(8) 


If  the  plant  described  by  Equations  (3)  and  (4)  contain^  integration 
operations,  they  nust  further  be  rt placed  by  numerical  integration 
approximations . 
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3.  Design  Problem 

The  original  PERSHING  weapon  system  attitude  control  system 
implemented  in  the  on-board  analogue  computer  may  be  represented  in 
block  diagram  form  (Figure  3),  where  <t> ,  <t> and  rj  represent  vehicle  pitc. 

attitude,  desired  pitch  attitude,  and  slant  altitude  velocity,  respec¬ 
tively.  Symbols  r  and  rj  represent  pitch  attitude  and  slant  altitude 

velocity  errors,  respectively,  where 


r(t)  =  $p(t)  -  <t»  (t) 


(9) 


and 


rje<t)  «  T)(t)  .  (10) 

In  this  report  only  the  pitch  attitude  control  system  will  be  redesigned. 
Specifically,  it  is  desired  to  match  the  continuous  system  signal  y(t) 
with  a  digital  control  system  y(t)  for  any  given  input  signal,  r(t).  Let 

y(t)  =  aQ  r(t)  +  x2(t)  ,  (11) 


whe^e 


x2(t)  =  £ 


a^sk's) 


(v  +  “)( 


’2s  + 


(12) 


Because  of  the  linear  characteristics  of  the  existing  continuous  system 
and  the  digital  system  to  be  proposed,  it  is  sufficient  to  match  the 
digital  system  variable  x2(t)  to  the  continuous  system  variable  x2(t). 

The  continuous  system  is  described  in  Figure  4,  where 


and 


d.  1 

a,  =  — 

1  T,T 


v  2 


a2  = 


d.  T1  +t2 


T 1T  2 


(13) 


(14) 
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In  the  state  space  formulation  of  Equations  (1)  and  (2), 


X  =  [Xj,  x2f  , 

r  .  [0,  rlT  , 

"o  1“ 

A  = 

-0  0. 

B  =  [0,  1]  , 

G(0)  =  [ai>  a2]  * 

E(0)  =  *,0^ 


(15) 

(16) 

(17) 

(18) 

(19) 

(20) 


The  equivalent  digital  system  wi:h  an  analogue  "plant"  (Figure  5)  is 
found  by  using  Equations  (5)  through  (8),  yielding 


G 


and 


where 


E (T)  =  a1a1  , 


and 


(21) 


(22) 


(23) 


(24) 
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If  the  integrators  of  the  analogue  plant  are  replaced  by  digital  numeri¬ 
cal  approximations,  D  ,  the  all-digital  control  system  is  of  the  form 

z 

shown  in  Figure  6  and  may  be  written  in  transfer  function  form, 


x2(z) 

R(z) 


E  (T)  D 


1  + 


(rw  \ 

“2  +  Dzal) 


(25) 


where  x^(z)  is  the  z-transform  of  the  output  x2(t)  of  the  all-digital 

system.  The  stability  considerations  brought  out  in  Reference  1  must 
be  adhered  to,  i.e.,  in  this  case  both  an  and  a2  should  be  positive 
numbers . 

3 

There  are  numerous  numerical  approximations  (D  )  of  integrators  . 

z 

As  might  be  expected,  the  simplest  formulations  give  the  poorest  results. 
Three  possible  numerical  integration  operations  are: 


! 


(Trapezoidal  Rule) 


(26) 


Dz  = 


-  t(i  +  4z_1  -r  z"2) 
-  — 


(Simpson's  1/3  Rule) 


(27) 


D  = 
z 


0.3  T  (l  +5z  ^-fz2  +  6z^+z^+5z^+z 


1  -  z 


-6 


(Weddle's  Rule) . 

(28) 


All  three  rules  have  essentially  ideal  phase  characteristics,  so  the 
quali:y  of  the  approximations  can  be  determined  primarily  by  considera¬ 
tion  of  their  amplitude  characteristics.  This  leads  to  the  ranking  of 


3 

Salzer,  J.  M.,  "Frequency  Analysis  of 
in  Real  Time,"  Froc .  IRE ,  Vol.  42,  February 


Digital  Computer  Operating 
1954,  p.  463. 
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Weddle's  rule  as  best,  Simpson's  1/3  rule  as  good,  and  the  trapez'ldal 
rule  as  poor.  However  the  trapezoidal  rule  is  simpler  to  program  and 

4 

introduces  attenuation  at  higher  frequencies. 


4.  Design  Problem  —  Numerical  Example 

Assume  that  representative  numerical  values  for  the  original 
PERSHING  control  system  are: 

aQ  =  -6,  aL  =  -2.4  .  (29) 

If  a  step  input  of  amplitude  R, 

r(t)  =  R,  t  2  0 

=  0,  t  <  0  ,  (30) 

is  applied  to  the  continuous  system,  the  x^(t)  response  is  formed  from 
Equation  (12): 


x2(t)  =£ 


-1 


W 


s  +  a2s  +  aj 


a-jC^Re 


2  2  sinj^  -  \a\  t 

A  -  \  i 


(31) 


Applying  the  final  value  theorem,  the  steady  state  value  x2(t)  is  zero. 

If  it  is  desired  to  obtain  values  of  f  and  u  of  0.707  and  84  radians/ 

J  n 

second,  resp^ativel’ ,  and  are  calculated  to  be 


ax  =  7056,  a2  =  118.8 


(32) 


L 

lou,  J. 
Book  Company, 


T.,  Digital  and  Sampled-Data  Control  Systems,  McGraw-Hill 
Inc.,  New  York,  1959  pp.  20f-207. 
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This  leads  to  a  time  response  of 


x2(t)  =  -287.4  Re"59,4tsin  58.92  t 


(33  ^ 


assuming  zero  initial  conditions,  which  is  plotted  in  Figure  7. 
Similarly,  if  a  ramp  input  of  slope  R, 

r(t)  =  Rt,  t  ^  0 


=0,  t  <  0 


(34) 


is  applied  to  the  continuous  system,  the  x^(t)  response  is 


x2(t)  =  £ 


-1 


a-o^R 


s^s2  +  a2s  + 


■  *iRLl  - 


-  {ajt 

e  sec  9  cos 


(/“ i  •  i  4  ‘ +  s). 


(35) 


^d.,.  -li 
0  =  tan 


a„ 


,  /  12 

2./al  -  4  aZ 


The  steady  state  v«-lue  of  x^(t)  is  seen  to  be 


Um  x  (t)  =  a,U 


t  — >  00 


(36) 


Again  using  the  representative  numerical  values  of  Equations  (29)  and 
(32),  the  ramp  r(  jponse  of  x2(t)  is  fc  ad  to  be 

x2(t)  =  -2.4  R[l  +  1.420  e”59-4tcos(58.92  t  -  0.7894)1  (37) 

which  is  plotted  on  Figur >  8  for  R  equal  to  unity. 
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The  proposed  digital  control  fystem  responses  to  step  and  ramp 
inputs  will  now  be  found  and  comp ^red  to  the  continuous  system  responses 
of  Equations  (33)  and  (37).  For  the  sake  of  computational  simplicity, 
Equation  (26)  (trapez  tidal  rule)  is  selected  fov  the  numerical  approxi¬ 
mation  of  integration.  Application  of  Equation  (26)  to  Equation  (25; 
leads  to  a  relation  for  obtaining  values  of  x2(nT)  at  sampling  instants: 


x2(nT)  =  Z 


|tf(t)( 

>R(z) 

(l  +  Ta2  +  ^  T2aL^z2  -  2 

(38) 


where 


R(z)  =  zfr(t)) 

If  the  final  value  theorem  is  applied  to  Equation  (38),  steady  state 
values  of  for  step  and  ramp  inputs  represented  by  Equations  (30)  and 

(34),  respectively,  are  found  to  be  identical  to  the  original  continuous 
system  steady  state  values  of  zero  and  a^R,  respectively. 

If  the  PERSHING  sampling  period  of  T  =  1/122  second  and  Equations 
(29)  and  (32)  are  applied  to  Equation  (38),  the  values  of  x2(nT)  may  be 
found  from 


x2(nT)  = 


■24.93j-r 


uT 


r(n  -  2)T.l  + 


1.314 x. 


-  0.4845  x. 


(n  -  1)T 


(n  -  2)T 


(39) 


”here  r(n  -  k)I 


and  x2 


(n  -  k)T 


denote  the  values  of  r(t)  and  x2> 


— 1  th 

x  pectively,  at  the  (n  -  k)  sampling  instant.  The  responses  to  unit 
step  and  unity  slope  ramp  inputs  are  calculated  from  Equation  (39)  and 
plotted  in  Figures  7  and  8  so  they  may  be  compared  directly  to  the  con 
tinvous  system  responses  to  the  same  inputs. 


To  show  the  effect  of  the  numerical  approximation  for  the  integra¬ 
tion  operation,  a  digital  formulation  of  the  control  system  using  ideal 
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integrators  is  analyzed  (Figure  5).  The  resulting  response  x2(nT)  that 
is  to  be  matched  to  the  original  x/>(t)  is  found  directly  from  Equations 
(3),  (4),  and  (21)  through  (24),  yielding 


x2(nT)  =  Z 


_ T(z  -  1)  E(T)  R(z) _ 

z.2  -  (l  -  a2T  -  |  ajT2  J  z  +  ^l  -  a2T  +  ~ 


(40 1 


Application  of  the  final  value  theorem  tqJJquation  (40)  yields  the  same 
steady  state  values  for  x2  as  for  x2  and  x2-  Again  using  T  =  1/122 

second  and  applying  Equations  (29)  and  (32)  to  Equation  (40)  yields  the 
equation 


V"T>  ■  -71-25[r(„  .  l)T  •  r(„  .  2)T]+  1-WuSj 

(41) 

The  responses  to  unit  step  and  unity  slope  ramp  inputs  are  calculated 
from  Equation  (41)  and  plotted  in  Figures  7  and  8  for  comparison  with 
continuous  system  and  all-digital  system  responses  to  the  same  inputs. 


0.3848  x. 


(n  -  1)T 


(n  -  2)T 


5.  Conclusions 

The  step  and  ramp  responses  of  the  original  continuous  analogue 
computer  control  system  and  the  proposed  all-digital  control  system  are 
shown  in  Figures  7  and  8  for  ease  of  comparison.  It  is  seen  that  the 
digital  system  response  x2  follows  the  continuous  system  response  r, 

with  a  time  lag  and  lower  peak  magnitudes.  The  digital  system  response 
can  be  made  to  better  match  the  continuous  system  response  if  a  more 
accurate  numerical  integration  routine  is  selected,  such  as  Equation 
J27)  or  (28).  The  besc  that  can  be  attained  is  indicated  by  the  response 
x^  of  the  digital  system  with  perfect  integrators  (Figures  7  and  8) ,  If 

a  better  match  of  responses  is  needed  usinji  the  technique  of  this  report-, 
the  sampling  period  will  have  to  be  decreased. 

The  technique  used  in  this  report  is  simple  and  easy  to  apply  and 
appears  to  yield  a  suitable  design  tool  for  designing  an  all-digital 
control  system  from  the  original  PERSHING  continuous  analogue  computer 
control  system.  This  statement  is  based  on  the  analytically  obtained 
results  presented  in  this  report.  It  is  strongly  recommended  that  the 
proposed  digital  implementation  be  checked  ou  on  a  computer  simulation 
of  the  PERSHING  attitude  control  system  to  confirm  the  results  reported 
here . 
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The  area  of  digital  redesign  of  a  continuous  control  system  appears 
to  be  fruitful  for  further  applied  research  and  poorly  understood.  To 
merely  add  a  digital  system  and  retain  the  same  numerical  gains  that 
were  used  for  the  continuous  system  will  not  provide  an  adequate  response. 


rftl 


m  (*) 


x  h)  ■  A  x  (t)  +  B  m  (t) 


x  (t) 


1 

G  (0) 

,  J 

Figure  1.  Continuous  Model 


Figure  2.  Sampled-Data  Model 


Figure  3.  Simplified  Block  Diagram  of  PERSHING  Analogue  Contiol  System 
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Figure  4.  Continuous  Pitch  Attitude  Control  System 


i.gure  6.  All-Digital  Model 


nse  (T  =  1/122) 


